arXiv:1506.00361vl [hep-lat] 1 Jun 2015 


The chromomagnetic operator on the lattice 


M. ConstantinouT M. Costa^, R. Frezzotti^, V. Lubicz^’^ 
G. MartinellT’T D. Meloni/ c,d ) H. Panagopoulos^, S. Simula^ 

Department of Physics, University of Cyprus, Nicosia, CY-1678, Cyprus 
Dip. di Fisica, Universitd di Roma “Tor Vergata” and INFN, Sezione di “Tor Vergata”, 
Via della Ricerca Scientifica 1, 1-00133 Rome, Italy 
Dip. di Fisica, Universitd Roma Tre, Via della Vasca Navale 84, 1-00146 Rome, Italy 
i d ) INFN, Sezione di Roma Tre, Via della Vasca Navale 84, 1-00146 Rome, Italy 
SISSA, Via Bonomea 265, 1-34136, Trieste, Italy 
(f) Dip. di Fisica, Universitd di Roma “La Sapienza” and INFN, Sezione di Roma, 

P.le Aldo Moro 2, 1-00185 Rome, Italy 



Abstract 

We present our study of the renormalization of the chromomagnetic operator, Ocm , which appears in the effective 
Hamiltonian describing AS = 1 transitions in and beyond the Standard Model. 

We have computed, perturbatively to one-loop, the relevant Green’s functions with two (quark-quark) and three 
(quark-quark-gluon) external fields, at nonzero quark masses, using both the lattice and dimensional regularizations. 
The perturbative computation on the lattice is carried out using the maximally twisted-mass action for the fermions, 
while for the gluons we employed the Symanzik improved gauge action for different sets of values of the Symanzik 
coefficients. We have identified all the operators which can possibly mix with Ocm , including lower dimensional and 
non gauge invariant operators, and we have calculated those elements of the mixing matrix which are relevant for the 
renormalization of Ocm- 

We have also performed numerical lattice calculations to determine non-perturbatively the mixings of the chro¬ 
momagnetic operator with lower dimensional operators, through proper renormalization conditions. For the first 
time the l/a 2 -divergent mixing of the chromomagnetic operator with the scalar density has been determined non- 
perturbatively with high precision. Moreover, the 1/a-divergent mixing with the pseudoscalar density, due to the 
breaking of parity within the twisted-mass regularization of QCD, has been calculated non-perturbatively and found 
to be smaller than its one-loop perturbative estimate. The QCD simulations have been carried out using the gauge 
configurations produced by the European Twisted Mass Collaboration with Nf =2 + 1 + 1 dynamical quarks, which 
include in the sea, besides two light mass degenerate quarks, also the strange and charm quarks with masses close to 
their physical values. 


1. INTRODUCTION 


A very natural explanation for the extraordinary success of the Standard Model (SM) in describing the electro- 
weak and strong interactions at the fundamental level is that the SM Lagrangian contains all relevant operators of 
dimension d < 4 composed by the (already observed) elementary particle fields and compatible with the principles of 
Lorentz invariance and gauge symmetry. The effects of higher dimension (d > 4) effective operators, which are not 
included in the SM Lagrangian, are expected to be naturally small, being suppressed by negative powers of the high 
energy scale M characterising the physics beyond the SM, as M 4 ~ d (up to logarithmic corrections). 

In this picture, a special role is played by the operators of dimension d = 5, as their contribution is suppressed 
by only one power of the high scale. In the leptonic sector, an important example of d = 5 operator is represented 
by the Weinberg operator for neutrino masses [l|, composed by two lepton doublets and two Higgs fields. After the 
occurrence of spontaneous electroweak symmetry breaking this operator leads to a natural explanation of the small 
light-neutrino masses, which are thus predicted to be inversely proportional to the large scale M. 

In the quark sector the d = 5 magnetic operators, which induce A F = 1 flavor changing transitions, are of relevant 
phenomenological interest. In the strangeness changing AS = 1 case, for instance, these magnetic operators contribute 
to both CP conserving and CP violating rare kaon decays, as well as to K° — K° mixing and to the CP violating ratio 
e'/e. In a large class of new physics models these contributions can be substantially larger than in the SM, which 
motivates the interest in studying their effects. This is the case, for instance, of generic supersymmetric extensions 
of the SM, in which AS = 1 transitions described by the magnetic operators are mediated by the strong interactions 
through virtual gluino exchanges 00 - 

For definiteness, let us consider in the following AS = 1 transitions. In both the SM and beyond, the low-energy 
effective Hamiltonian contains four magnetic operators of dimension 5, 

h as=i, d =5 = + ciQi) + h.c. (i) 

2 =± 

which are defined as: 

($sL <r» v dR ± $ sR a^ ip dL ) , ( 2 ) 

Qf = G ^^dR±i’sRCT l * V ■ (3) 

In the above expressions, F^ v and G^ v represent the electromagnetic and strong field strength tensors respectively, 
ip s and V'd are the strange and down quark fields and the subscripts R , L denote the left/right chiral structure (l±q 5 ). 
The coefficients (7* and C®, multiplying the electromagnetic (EMO) and chromomagnetic (CMO) operators in the 
effective Hamiltonian, contain the effects of the physics at short distance and they depend on the specific structure 
of the new physics model. These coefficients can be calculated perturbatively via the OPE. The long distance effects 
of the strong interactions are encoded in the operator matrix elements and thus require, for their evaluation, a 
non-perturbative method, primarily a lattice QCD calculation. 

The matrix element of the EMO between kaon and pion states may be relevant in the CP-violating part of the 
K —> 7 t£ + £~ semileptonic decays (see Ref. 1) and its determination offers, for instance, the possibility to put bounds 
on the supersymmetric couplings related to the splitting of the off-diagonal entries in the down-type squark mass 
matrix. The matrix element (7r|Q+|IF) has been computed on the lattice both in the quenched approximation Q and 
with Nf = 2 flavors of degenerate sea quarks [g]. 

Several matrix elements of the CMO between kaon and pion states are of phenomenological interest for supersym¬ 
metric extensions of the SM. The matrix element (7r 0 |Q+|AT 0 ) may provide contributions to the K° — K° mixing 
amplitude (see Ref. 0), while the matrix element (7r + 7r | Q g \K°) may play a role for determining s '/e and for the 
A I = 1/2 rule (see Ref. 0 ). Finally the matrix element (7r + 7r + 7r _ |Q+| K + ) may contribute to the CP-violating part 
of the K —> tttttt decays [3(. All the above-mentioned matrix elements can be parameterized in terms of suitably 
defined B-parameters: 


(AQj\ K °) 

(n+n-\Q-\K°) 

(7r + 7r + n-\Qj\K+) 


_L m k (PK ' Ptt) 

v^2 327r 2 m s + rrid CM 

■ 11 Mfc Ml „K2 tt 

32tt 2 (m s + m d ) CMO ’ 

11 Ml Ml nK3n 

167r 2 fl (m a + m d ) CMO ' 


( 4 ) 


At leading order (LO) in Chiral Perturbation Theory (ChPT) the CMO has a single representation in terms of 
pseudo-Goldstone boson fields 0. Therefore, the three B-parameters appearing in Eq. © are related by chiral 
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symmetry, which predicts at LO their equality: 

Bcmo = Bcmo = Bcmo — Bcmo ■ ( 5 ) 

A lattice calculation of the matrix elements of the CMO is challenging, particularly when more than one pion are 
considered in the final state. Even in the case of only one final pion, which corresponds to the matrix element 1 K) 

of the operator Q+ = (g/ 16n 2 ) ip s a^ G^if^d, n o results have been produced so far. The main difficulty, with respect 
to the EMO case, is that strong interactions induce a mixing of the CMO with operators of lower dimension, with 
coefficients which are power divergent with the cutoff, which on the lattice is the inverse of the lattice spacing 1 /a. 
The leading divergence of the bare CMO, which is of order 1/a 2 , is determined by the mixing with the dimension-3 
scalar operator V’sV'd- Its coefficient must be evaluated in a fully non-perturbative way, since non-perturbative effects, 
e.g., factors of the form aA qcd, combined with factors which diverge as inverse powers of the lattice spacing can give 
finite (or even divergent) contributions §. 

In order to define the properly renormalized CMO, besides the subtraction of the lower dimension operators, the 
mixing with equal dimension (d = 5) operators, including the CMO itself, must also be taken into account. This 
mixing is only logarithmically divergent and can be thus evaluated in perturbation theory. The one-loop calculation 
of the corresponding renormalization factor and mixing coefficients is one aim of the present study. Specifically, we 
have considered a lattice regularization of QCD defined by a generic class of Symanzik improved gluon actions and 
a twisted-mass quark action ins. By investigating the symmetry properties of this action, we have shown that 
the renormalized CMO mixes with a total of 13 operators (including itself), of which seven are not present on-shell; 
among them, there will be non-gauge invariant (but BRST invariant) operators as well. For on-shcll matrix elements, 
the mixing assumes the general form: 

_ _ 

goi’s & [AIS Gfj.ui’d = Zi [g 0 ip s & [IIS G^ d ] + Z 2 [{m 2 d + m «)V’«V , d] + z 3 [ m d m s il) a il) d ] 

_ ££ _ 

+Z 4 [□ (^ d )\ + ^12 [ % (j'd TH'd H” ^s s'ys'^d] + Z i3 [V> s ip d \ , (6) 

where R denotes the corresponding renormalized operators, Z 42 oc 1/a, Z 13 oc 1/a 2 , r s t d \ is the Wilson hopping 
parameter of the strange (down) quark, and we have evaluated the renormalization factor Z\ and the mixing coefficients 
Z 2 — Z\z at one loop. Note that the presence of the mixing with the pseudoscalar density is due to the parity violation 
in the twisted-mass formulation of QCD on the lattice. As discussed above, the power-divergent coefficients Z\ 2 and 
Z 13 require an independent non-perturbative determination. 

A strategy to implement non-perturbatively the subtraction of the mixings of the CMO with the lower dimension 
operators has been anticipated in Refs. 1 1H13| . In this work we apply it for obtaining the first non-perturbative 
determinations of the power-divergent mixings Z13 and Z\ 2 . To this end we have used the gauge configurations 
produced at three values of the lattice spacing (between ~ 0.6 and ~ 0.9 fm) by the European Twisted Mass 
Collaboration (ETMC) with Nf = 2 + 1 + 1 dynamical quarks, which include in the sea, besides two light mass 
degenerate quarks, also the strange and charm quarks with masses close to their physical values EMU. It turns 
out that the one-loop perturbative estimate of Z13 differs only by less than 10 % from the non-perturbative results 
at the three values of the lattice spacing, while the non-perturbative determination of Z 42 is found to be smaller 
than the corresponding one-loop perturbative result. This finding suggests that, together with our non-perturbative 
determinations of the mixing coefficients Z 13 and Z 12 , the perturbative estimates of the renormalization factor Z\ 
and of the mixing coefficients Z 2l Z3 and Z 4 may be used for the determination of the (renormalized) CMO matrix 
element. Preliminary results for such a matrix element between pion and kaon states have been presented in Ref. G2 
and the final ones will be the subject of a forthcoming publication [l 8 j . 

The outline of this paper is as follows. Section [2] provides a brief theoretical background in which we introduce 
the symmetries of the employed actions and the transformation properties of all candidate operators which can mix 
with Ocm at the quantum level. Section [3] contains a summary of the computational procedure for the Green’s 
functions of the chromomagnetic operator. This Section is divided in two subsections. In Subsection 13.II calculating 
the 2-point and 3-point Green’s function of Ocm in dimensional regularization (DR), we construct a set of eleven 
independent equations for the disentanglement of the mixing coefficients. We present these coefficients in the MS 
renormalization scheme. On the other hand in Subsection 13.21 using the lattice formulation and the results which we 
found in Subsection 13+1 we calculate the mixing coefficients on the lattice, again in the MS renormalization scheme. 
In Section [|] we describe the first non-perturbative, high-precision determination of the 1 /a 2 -divergent mixing of the 
chromomagnetic operator with the scalar density, using the ETMC gauge configurations with Nf = 2 + 1 + 1 produced 
at three values of the lattice spacing. We also describe the first non-perturbative calculation of the 1/a-divergent 
mixing of the chromomagnetic operator with the pseudoscalar density using the renormalization condition introduced 
in Ref. (Tlj| . Finally, we conclude in Section [5] with a discussion of our results and possible future extensions of our 
work. For completeness, we have included two Appendices containing the mixing coefficients Z.- t (Appendix [A| and 
the one-loop perturbative renormalization factors Z c , Z^, Z mi Za and Z q on the lattice (Appendix iBl). Preliminary 
results for the above coefficients have been already presented in Ref. [1 1! Il.'ll . ft7l ] . 
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2. SYMMETRIES OF THE ACTION AND TRANSFORMATION PROPERTIES OF OPERATORS 


We start by studying the mixing of the chromomagnetic operator 1 : 

Ocm = go 4> s CT/ii/ G^d, (7) 

using both DR and lattice regularization (L). On the lattice we use the fermion setup studied in Refs. urn ; in 
particular, valence quarks are described by the twisted-mass/Osterwalder-Seiler (OS) action at maximal twist, which 
in the physical basis reads: 


where 


S F [^f,^f,U] = « 4 EE 1>f (x) 

f * 


7 • V - *75 W cr (r/) + rrif ipf(x) 


7 ‘ V = i^CV^ + V^), 
w cr (r f ) = ^E v ;v,+M cr ( r/ ) : 


( 8 ) 


(9) 

( 10 ) 


r/ is the Wilson parameter for the flavor f = u, d, s and Af cr (r/) is the corresponding critical quark mass. 

The full fermion action includes also a part describing sea quarks, and possibly a ghost part (to compensate the 
valence quark determinant for the partially quenched flavors) [l9(; these parts will not be needed in our perturbative 
calculation. For the gluon part we employ the Symanzik improved action: 




^ [ c 0 J2 Re Tr {1 “ tfpiaq.} + ci ^ Re Tr {1 f^rect.} 

0 plaq. rect. 

+ C 2 J2 Tr { 1 - t/chair} + C 3 J2 Tr { 1 ^paral.}j 5 

chair paral. 


(ii) 


where the Wilson loops are products of consecutive links in the directions (p,u,—p,—u), (p,v,u,—p,—u,—u), 
(p,u,-p, p,-i/,-p), {p,v,p,—p,—v,—p) for t/piaq.,d r ect.,t4hair and C/parai., respectively. The Symanzik coefficients 
cq,ci,C 2 ,C 3 may take arbitrary values, subject to the constraint: 


c 0 + 8 ci + I 6 C 2 + 8 C 3 = 1 , 


( 12 ) 


which ensures the correct classical continuum limit. Our results (Section 13.21) will be provided for some of the most 
popular choices for Cj. 

There exist certain symmetries of the action (valid both in the continuum and lattice formulation of the theory) 
which reduce considerably the number of operators that can possibly mix with Ocm at the quantum level. These 
symmetries are defined by means of the discrete transformations V (continuum parity), 

( Uo(x) -A U 0 (x v ), Uk{x) -A Ul(x-p — ak ), fc = 1,2,3 
V : < ifjf (x) —> 70 tpf{xv) (13) 

{ -t ipfix-p) 7o , 


where x-p = (—x, Xq) and p is the unit vector in the //-direction, 


f Up(x) 

V d : l f f {x) 


u l{- x ~afi) 

e 3 z 7 r / 2 i(jf(—x) 

e 3 ™/ 2 ij) f (—x), 


{ i>fi x ) 

which, besides inverting x —> —x, counts the parity of the dimension d of each operator by multiplying it by e™ d , 


(14) 


77-5 = 17 77./; 
/ 


77/5 : 


75^/ 

~7/75 , 


(15) 


1 In our notation go is the bare coupling constant, il) s ,d are the s- and d-quark fields, is the gluon tensor and cr^ v = (i/2) [ 7 ^, 7 ^]. 
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C (charge conjugation; T means transpose) 

ip(x) -*■ -if>{x) T i'y 0 j 2 (16) 

U^x) U*{x ), m = 0,1,2,3, 

and 5 (exchange between the s and the d quark) 

{ Vb(z) •O- ipd{x) 

4> s {x) O ^ d (x) (17) 

m s -O m d . 

In terms of the above transformations, the symmetries of the action are 2 : 

• V x T>d x (m —^ —m), where m are all masses except M cr 
• Dd x 

• C x S, if r s = r<j (18) 

• C xV x S, if r s = —ru . 


In order to identify which operators can possibly mix with Ocm , we examine the transformation properties of all 
candidate operators under the above symmetries; admissible operators must transform in the same way as Ocm- 
Furthermore, by general renormalization theorems, these operators must be gauge invariant, or else they must vanish 
by the equations of motion. 

In TableUwe present all candidate operators along with their transformation properties. Operators marked by " 
have the same properties as Ocm and thus may mix with it. Operators marked by "(+)" or "(—)" have the same 
transformation properties as Ocm only if r s = r d or r a = —r d , respectively; for this reason the Wilson parameters 
r a , r d have been explicitly introduced in On and On below (see Eqs. ® - ®)- There follows immediately that 
Ocm = 0 i can only mix with the following operators: 


01 

= goiPsO^G^ipd 

(19) 

02 

= (rn 2 d + m 2 s )tp s ilj d 

(20) 

03 

= m d m s ip s i/jd 

(21) 

04 

= 

(22) 

05 

= ip s (-^P + m s )(0 + m d )4>d 

(23) 

06 

= %Ij s C 0 + m d ) 2 ip d + + m s ) 2 ip d 

(24) 

07 

= m a ij; s (0 + m d )ip d + m d ip s CP + m s )ipd 

(25) 

08 

= m d ip s ($ + m d )ilJd + rn s il> s C < p + m s )4> d 

(26) 

09 

= ip s ^0 +m d )^d ~ i> a {-*P + m 8 )$ipd 

(27) 

010 

= ip a ^(0 + m d )i> d - tl> a {-P + m s )]!)ip d 

(28) 

011 

= ir d ip a j 5 (-0 + m d )4>d + ir a i+ rn^^ipd 

(29) 

012 

= i(.r d m d +r s m s )'il> s 'y 5 ip d 

(30) 

013 

= 

(31) 


where □ = and the left and right covariant derivatives are defined in terms of the gluon field as follows: 

n = & n T igoA^ , ^ = d [i ig^A^ . (32) 


2 Note that, in the case of r s = — CVS will not be a symmetry of the valence part of the action which contains a u quark, since it 
will require r u —> —r u - However, the u quark can be dropped from the valence part of the action, since our operator does not contain u 
quarks, and therefore the Green’s functions of interest will also not contain any external u quarks. Nonetheless, it is important to note 
that the sea quark part of the action is symmetric even in the presence of u, since it is an even function of the Wilson coefficients rf 
(by virtue also of M C r(—rf) = —M C r(rf)) 0. 
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Operators 

V x T> d x 

'Dd x 7^-5 

CxS 

CxVxS 


{m —> —m) 


if r a = r d 

if r s = -r d 


Dimension 3 operators 


/ 


- 

+ 

+ 

+ 


iip^ipd 

+ 

+ 

+ 

- 


Dimension 4 operators 



(m d + m s )ip a tl> d 

+ 

+ 

+ 

+ 


{m d - m a )ip 3 il) d 

+ 

+ 

- 

- 

(+) 

i(m d +m s )ip 3 'y5ip d 

- 

+ 

+ 

- 

(-) 

i (m d - ra s )i/i,7s^d 

- 

+ 

- 

+ 


+ m d )ip d + + m s )ip d 

+ 

+ 

+ 

+ 


ipstfi + m d )ip d - tp s {-'p + m a )ip d 

+ 

+ 

- 

- 

(+) 

+ m d )tp d +iif> s (—'p + m a )-y 5 ipd 

- 

+ 

+ 

- 

(-) 

+ m d )ip d - iip s (—'p + ms^sipd 

- 

+ 

- 

+ 


Dimension 5 operators 


/ 

go ip s <j iii/GLii/ipd 

- 

+ 

+ 

+ 


i go ips'ysGfiuGd 

+ 

+ 

+ 

- 

/ 

{m 2 d + m 2 3 )ip s ip d 

- 

+ 

+ 

+ 


i {ml + m 2 a )'i/j s 'y 5 'il! d 

+ 

+ 

+ 

- 


{m 2 d - ml)if> s il) d 

- 

+ 

- 

- 


i {m 2 d - m 2 s )i/j s -y 5 ilj d 

+ 

+ 

- 

+ 

/ 

m d m 3 ip s ipd 

- 

+ 

+ 

+ 


im d m s tp g ^ 5 tp d 

+ 

+ 

+ 

- 

/ 

m a ip a (P + m d )ipd + m d 4’ a {-J) + m a )ip d 

- 

+ 

+ 

+ 

/ 

m d ip s (j5 + m d )ip d + m a 'ip a {-J) + m a )ip d 

- 

+ 

+ 

+ 


m a ii a (p + m d )ip d - m d ip s {-'p + m a )ip d 

- 

+ 

- 

- 


m d ip a (j5 + m d )ip d - m a ip 3 {-'p + m a )ip d 

- 

+ 

- 

- 


imsip^nCfi + md)ip d + im d ip a {-'p + m^sipd 

+ 

+ 

+ 

- 


im d ip s 75(^5 + m d )ip d + im a tp s {-'p + m a )j 5 ipd 

+ 

+ 

+ 

- 


imsip^sCfi + m d )ip d - im d ip a {-J) + m a )^5tpd 

+ 

+ 

- 

+ 


imdip^sCfi + m d )ipd - im a il) 3 {-J) + m a )j 5 il>d 

+ 

+ 

- 

+ 

/ 

ip a Cfi + m d ) 2 ip d + i} a {-Jj + m a ) 2 ip d 

- 

+ 

+ 

+ 


ip 3 Cfi + m d ) 2 ip d - ip a {-7) + m a ) 2 ip d 

- 

+ 

- 

- 


i ’'PslzCfi + m d ) 2 ip d +iip a {-J) +m a ) 2 'y s ip d 

+ 

+ 

+ 

- 


iip 3 -fs{J5 + m d ) 2 il) d ~ iip 3 {-'p + m a ) 2 -f 5 tpd 

+ 

+ 

- 

+ 

/ 

□ (V’sV’d) 

- 

+ 

+ 

+ 



+ 

+ 

+ 

- 

/ 

ip s {~'P + ms) (ft + m d )ip d 

- 

+ 

+ 

+ 


iip 3 {-J) + m a ) 75 C0 + md)ipd 

+ 

+ 

+ 

- 

/ 

+ rn d )ip d - ip 3 {~'P + m s )~$ip d 

- 

+ 

+ 

+ 

/ 

tpjtCP + m d )ipd - ip s {~7> + msfpipd 

- 

+ 

+ 

+ 


•tpJUfjS + m d )%pd + ij) 3 {-J) + m s )~^Tp d 

- 

+ 

- 

- 


+ m d )ip d + tpsi-'P + m s ))9ipd 

- 

+ 

- 

- 


iipj® jsCfi + md)ipd -iip 3 {-'P + m s )j 5 ~fiip d 

+ 

+ 

+ 

- 


+ rn d )ip d -iip s {-J) + m^-y^ip d 

+ 

+ 

+ 

- 


iipj&'ysffi + m d )ip d + iip 3 {-'p + m a ) 75 '^^ 

+ 

+ 

- 

+ 


i ipjl'ysCfi + m d )ip d + itp a {-^P + ipd 

+ 

+ 

- 

+ 


TABLE I: Transformation properties of gauge invariant opera|prs and of operators which vanish by the equations of motion, 
in the physical basis. 
































































We do not impose in our calculation the conservation of external momentum. Therefore, the list of independent 
operators in Table |T| accounts for operators which are total derivatives 3 . As for the parameters r s and r^, in our 
perturbative calculation we make the (independent) choices of values r s = ± 1 , = ± 1 , consistently with their values 

in the simulations. 

Operators O 9 and O 10 are not gauge invariant, but they are admissible candidates for mixing, since they vanish by 
the equations of motion; indeed, they will mix with Ocm both in DR and on the lattice. The operators On, 0 12 , 0 13 
are of lower dimension and thus they do not mix with 0 1 in DR; they do however show up in the lattice formulation. 

Before closing this Section we mention that in the presence of the electromagnetic interactions the operator Ocm 
can mix also with the EMO (see Eq. ©). The corresponding mixing coefficient is of order 0{g 2 ), i.e. it does not 
vanish formally in the limit of zero quark electric charge, since the latter is already included in the definition of the 
EMO. 


3. RENORMALIZATION FUNCTIONS 

The operators O R are related to the bare ones, Oi (i = 1,..., 13), through: 

13 

Oi = ^ ZijOf (in matrix notation : O = ZO R ), (33) 

1=1 

x y 

where the 13 x 13 mixing matrix Z, 3 (which should more properly be denoted as Z-' , where X = DR, L, ... is the 
regularization and Y = MS, RI' ,... the renormalization scheme) obeys: 

Z=t + 0(g 2 ), (34) 

where g is the renormalized coupling constant. Since we are interested in O r we only need to calculate the first row 
of the inverse mixing matrix, Z~ y , which, to one loop, is immediately related to the first row of Z: Zi = Zn. 

Since renormalization conditions are typically imposed on amputated renormalized Green’s functions, let us relate 
the latter to the bare ones. For the quark-quark Green’s functions: 

13 

= (z ^ 1 

i =1 
13 

= amp, i> = y/Z^1p R . (35) 

i=l 

The one-loop Feynman diagrams contributing to ('i/jO i V^amp are shown in Fig. [T] Note that Eq. ([35]) holds for an 
arbitrary regularization and arbitrary renormalization scheme; the only condition on the renormalization scheme is 
that it be mass-independent, in which case the quark field renormalization constant Z^ does not depend on flavor. 
To avoid heavy notation we have omitted coordinate/momentum arguments on O , as well as Dirac/flavor indices 
on (?/> O 'll )), etc. 


FIG. 1: One-loop Feynman diagrams contributing to the 2-pt Green’s function of the chromomagnetic operator, Oi . A wavy 
(solid) line represents gluons (quarks). A crossed circle denotes the insertion of G \. 


3 Instead of the operator O 4 one can consider the operator 0' 4 = \}j S *D which shares the same transformation properties of O 4 . It 
can easily be shown that the operator 0' 4 is a linear combination of O 4 and other operators entering the basis (1191311 ) . The choice of O 4 
has the advantage that it does not contribute to physical amplitudes, since it is a total four-derivative. Moreover, its non-perturbative 
determination is quite simpler, because its matrix element between physical hadron states is simply given by the corresponding matrix 
element of the scalar density multiplied by the squared four-momentum transfer. Finally, as it will be shown later, its mixing 

with the chromomagnetic operator is vanishing at one-loop and therefore any redefinition of O 4 does not change the mixing of the other 
operators of the basis at one-loop. 
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Similarly for quark-quark-gluon Green’s functions we have: 

13 

(V; ii Of^Af) amp = ^4 /2 ^(Z- 1 ) li (V>O^A !/ ) amp , A v = Jz~ a A r . (36) 

2=1 

Strictly speaking, on the right-hand sides of Eqs. (l35|) and (l36l) one must take the regulator to its limit value (i.e. 
e —> 0 in DR or a —> 0 on the lattice). This limit is convergent, provided all renormalization functions Z have 
been appropriately chosen. It is only in this limit that the right-hand sides of Eqs. © and ( 1 M 1 ) are equal to the 
corresponding left-hand sides. 

The one-loop Feynman diagrams contributing to {ip 0\ ip ^4„) amp are shown in Fig. [2] (one-particle irreducible (1PI)) 
and Fig. [3] (one-particle reducible (1PR)). While 1PR diagrams do contribute to the renormalized Green’s functions 

{ip R 0 R ip A R ) amp: their contribution to the mixing matrix cancels out. 



FIG. 2: 1PI Feynman diagrams which contribute to the 3-pt Green’s function of 0\. Diagrams 1, 4, 6 do not appear in DR. 
Wavy (solid) lines represent gluons (quarks). Crossed circles denote the insertions of (D\ . 


yzy 1 (~\ 

FIG. 3: 1PR Feynman diagrams which contribute to the 3-pt Green’s function of 0\. Wavy (solid) lines represent gluons 
(quarks). Crossed circles denote the insertions of 0\ . 

Imposing renormalization conditions of the above 2- and 3-pt Green’s functions is sufficient 4 in order to obtain all 

Zi- 

In some definitions of Ocm (see, e.g., | 20 |) there is an extra factor of a quark mass: 

Ocm = m Ocm , (37) 


4 One could of course calculate also 4-pt Green’s functions; in doing so, a number of consistency checks would emerge regarding the 
divergent part of the mixing coefficients Zi. Further Green’s functions (5-pt and above) will bring in no superficial divergences, and no 
further renormalization conditions (or consistency checks) will arise. 

































where m is the mass of one of the quark flavors. The renormalized mass m R is given by m R = Z^ m; in a mass- 
independent scheme, Z m is also flavor independent, like Z^. In this case the renormalization matrix Ztj for Ocm is 
simply given by: Z ig = Z m Zij. 

By analogy with Z m , a multiplicative factor of Z g must be included in Z-\ , if the calculation of Green’s functions 
involves the operator 4’ s a^ v G llu 'ilj ( i, rather than g^gCr^G^ipd- We will make use of this fact in Eq. (l58l) . The 
calculation of Z m and Z g is presented in Appendix [B] 

In order to impose renormalization conditions, we need the expressions for the tree-level 2-pt and 3-pt 
Green’s functions of Oi, i = 1,..., 13. The tree-level parts of the 3-pt amputated bare Green’s functions 
{ip s (q 2 )Oi(x)ijj d (q 3 ) A I/ (( 7 i)) amp are shown (apart from an overall factor of e lx 'G q i- q 2+ q 3)) j n Table HD similarly for 
the tree-level parts of the 2-pt bare Green’s functions {4’s( < l2)Oi(x)'i/j d (q3)) anip . Note that the tree-level 3-pt Green’s 
functions, despite being amputated, receive also contributions which are not 1PI, as shown in Fig. [4j We do not 
include these in Table El however, their value can be easily deduced from the corresponding tree-level 2-pt Green’s 
functions. 


FIG. 4: 1PR Feynman diagrams contributing to the tree-level 3-pt Green’s functions. Wavy (solid) lines represent gluons 
(quarks). Crosses denote the insertions of the operator Oi , i = 1,..., 13. 


3.1. Dimensional Regularization 

The next step in our renormalization procedure is to calculate the MS-renormalized 2-pt and 3-pt Green’s functions 
of Ocm ; in order to do so, we must regularize the theory in D-dimensions (D = 4 —2 e), in the continuum. The general 
form of the O (1/e) part of the bare Green’s functions is (consistently with the tree-level values of the operators in 

Table HD: 


\DR 


(V’C’l V’) amp 


1/e 


Pi (92 + 93 ) + P 2 ( m 2 s + m 2 d ) + p 3 i ( m d + m s 92 ) 


+P4 i (m s t/ 3 + m d 92 ) + P5 92-93 + p6 92 93 + P7 m s m d 


(38) 


(i/jO^A^ip! 


l/e 


Ri g (92 + 93)^ + R2 9 ilu & + 9^2 7v) + R3 i g ( m s + m d ) + R4 (-2 * g a pi/ q lp ) ( 39 ) 


where g is the renormalized coupling constant in the MS scheme, which is related to the bare coupling constant in 
DR, 9o R , through: g = p~ e (Z^AMS)- 1 g D R anc j Ru jj. are numerical coefficients. Additional terms in Eqs. (13811 
and (1391) (such as gqi u or g/h^ v in Eq. (15^1) ) would imply mixing with further operators from Table Q] but this is 
excluded by the symmetries listed in Eq. (Ell- 
Computing pi, Ri to one loop we find: 


Pi 


P 2 


P3 


9 


7 r 2 e 
2 1 


(-3 C F ) 


16 7 T 2 e 

q 2 1 


- (~6C f ) 


16 7T 2 


(3 C F ) 


Pi = P5 = P6 = P7 = 0 


(40) 

(41) 

(42) 

(43) 


9 










Operators 

Tree Level 2-pt 

Tree Level 3-pt (1PI) 

Or 


0 

2igoa 

02 

(m 2 d + m 2 s )ip s ip d 

m d + m 2 s 

0 

03 

m d m s ip a il)d 

rridm s 

0 

04 

□ (V’sV’d) 

(<72 — Q3) 2 

0 

05 

tp s {-'p + m a )(J) + md)tpd 

-0293 + i&m d + *9377*3 + m a m d 

-90(^27^ + 7*^93) + igo(m s + m d )-yu 

06 

ip s Cfi + m d ) 2 ip d + + m 3 ) 2 ipd 

-92 - 93 + 2i(md93 + m s q 2) +m 2 d + m 2 s 

—2(/ 0 *(7 M ^qi M - 2go(93i- + 92*/) — 2ig 0 (m d + 77*3)7!/ 

07 

m a ip s (j5 + m d )^d + rrid'ipsi-'fi + m s )ipd 

m a ( *93 + m d ) + m d (iq 2 + m s ) 

igo(m s + md)7„ 

08 

m d ip s (f!) + m d )ipd + m s ilj a (-'p + m a )ip d 

m d { *93 + m d ) + m s (iq 2 + m a ) 

igo(m s + m d ) 7*/ 

09 

+ md)^ - V> s (-^ + m a )~Jfif> d 

~ *(92***d + 93«*s) 

90(9271/ + 

010 

+ md^d - ipsi-'P + ms^tpd 

-93 - 92 + i{(/ 2 m a + 93 m d ) 

-2ga{q2 V + 93!/) + 9o(9Wi/ + 7^93) - 2*g 0 <7 M i/9i*i 

011 

+ md)tl’d + ir 3 ip s {-'p + m a )^d 

ir d j5 (*93 + ***d) + *77 (*92 + m a ) 75 

-go(r d - 77)757*/ 

012 

i ( r d md + r s m s )ip a 'y5ipd 

* (*•<* m d + 77 77*3)75 

0 

013 

ip s ipd 

1 

0 


TABLE II: Operators which will possibly mix with the chromomagnetic operator in the physical basis, along with their tree-level 2-pt and 3-pt (1PI) Green’s functions. Here, 
qi is the external gluon momentum and 52,3 is the external final (initial) quark momentum. 


















(44) 


Ri = 
R 2 = 

R 3 = 

i?4 = 


9 


1 


(—6 C F ) 


16 7r 2 e 

9 2 1 {3 Nr. 

16 7t 2 e 
■ 2 1 


9 


16 7T 2 e 
■ 2 1 


3 N r 


2 N r 


9 


1 


16 7 T 2 e \ Nr 2 Nr 


7 N c 3 aN c 

- T~ + —^ 


(45) 

(46) 

(47) 


Here, N c is the number of colors, Cf = (-/V 2 — 1)/(2 7V C ) is the quadratic Casimir operator in the fundamental 
representation, a is the gauge parameter (a = 1(0) corresponds to Feynman (Landau) gauge). 

We have also computed the finite parts (O(e 0 )) for the above Green’s functions, which are just the corresponding 
MS-renormalized Green’s functions. These are irrelevant for the computation of the mixing coefficients in the MS 
scheme in DR; however, they are necessary in the calculation of Z r3 with lattice regularization and MS renormalization, 
see Section lT2l Using the form of Eqs. (l38l) - (l39ll and the tree-level Green’s functions of the various operators (Table 
ED), we construct a set of equations for the disentanglement of the mixing coefficients; in particular, by demanding 
that the coefficients of 0( 1/e) in the left-hand sides of Eqs. (13511 - (l36l) vanish 5 , we obtain, to order g 2 : 


7 dr,ms 


ryDR, MS ryDR, MS 

^6 Z 10 

= Pi 

(48) 

ryrDRMS . ryDR, MS . ryDR, MS 
^2 + + ^8 

= P2 

(49) 

r, ryDR, MS . ^DR, MS . ryDRMS 

Z z-g + Zj s + ^10 

= P3 

(50) 

ryDR. MS . ryDRMS 

Z,g Zj 7 + Zj g 

= P4 

(51) 

ryDR,MS 
— Zr 4 

= P5 

(52) 

ryDR, MS 2 r^DRMS 

= P6 

(53) 

ryDRMS ryDRMS r, ryDRMS 

— Zvg — Zvg — Z Zrj 

= P7 

(54) 

7 DR,MS r, ryDR, MS r, ryDR, MS 

Z Z Zj iq 

= Ri 

(55) 

ryDRMS . ryDR,MS . ryDR,MS 
Z* 5 + ^ 9 + ^10 

= r 2 

(56) 

r, ryDR, MS . ryDRMS . ryDRMS 
— z Zj g + Zj 7 + Zr s 

= r 3 

(57) 

2 DR. MS . ryDR,MS . ryDRMS 

g z 1 +z 6 ’ +z 10 

= Ri+g 1 +o Z A +Z 9 j 

(58) 


where zf >R,MS = 1 -| - g 2 Q(g 4 ) and -2q>i = 0(g 2 ) by Eq. (1M1) : also: 


ryDR, MS -i | Z 

z \b =l + g z. 


7 dr ,ms _ 


= 1 + 9 s z_ 


DR, MS 

+ 0(0), 

DR, MS 

1 

1 

'0 

z i> 

16 7 r 2 

e 

DR, MS 

+ 0(g 4 ), 

DR,MS 

1 

1 

A 

Z A 

16 7T 2 

e 

,DR, MS 

+ 0(g 4 ), 

DR,MS 

1 

1 

'a 

z g 

16 7T 2 

e 


iN c 2 N f 


6 2 

U 11 N r . 


(59) 

(60) 
(61) 


In particular, Eq. 


3 6 

stems from the requirement that the coefficients of [—2iga^ v q\ ii )lt in the left-hand side and 


5 Note that Eq. (1361) will also contain 0( 1/e) terms which are not polynomial in qi , m; such terms arise from the 1PR one-loop 3-pt 
Green’s function of 0\ (Fig. [3J and from the 1PR tree-level Green’s functions of O 2 ,..., 0 13 (Fig. QJ. By Eq. (1.351) all such terms cancel 
out among themselves. 
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right-hand side of Eq. (13^1) coincide: 


0 = (1 + g 2 z, 


2 DR, MS\ 


MKn /b 1 o DR MS 




1 + - 2 Z DR,MS 

2 


)( 


1 _L „2 DR,MS 

1 + g Z 


) ^ 


2 D-R,MS\/-> i d \ r 7 
- g z i )(i + i? 4 ) -z< 


Dfl.MS 


- z 


Dfl.MS 

10 


(62) 


only the 0( 1/e) part 


As it stands, the system of 11 equations (Eq. (1151) - (155} 1 for the 10 unknowns Z® R,MS — Z R a R ' MS appears over¬ 
constrained; indeed, Eqs. (l48l) . (1521) and (l55l) can only be compatible if 2 p\ = R\. This relation is indeed confirmed 
by our results (Eq. (l40ll and Eq. (jj4j| ). The presence of Z RR ’ MS in Eq. (l58l) stems from the fact that all one-loop 
Green’s functions were calculated with an insertion of ip s a^ Gipd (rather than g ip s G^ ipd ; thus multiplication 
by Z g is necessary in a way analogous to Eq. 437}. 

Solving the above equations, we obtain the mixing coefficients: 


ryDR. MS 


ryDR, MS 
Z 2 


ryDR ,MS 
Z 3 

ryDR, MS 


Z 5 


ryDR, MS 
Z 6 


ryDR, MS 

z 7 


ryDR, MS 
Z 8 


ryDR, MS 
Z 9 


ryDR ,MS 
^10 


1 + 


16 7T 2 


1 

e 




9 


16 7T 2 

0 


- iN ' + K 


0 

g 2 1 /3 Nc 3 \ 

16 7 r 2 e \ 2 N c ) 

0 

g 2 1 / 3JV e 3 \ 
16 7 T 2 e V 4 + 2 No) 


0 

9 2 1 (3Nc _ 3 \ 

16 7 T 2 e \ 4 2 N c J 

g 2 1 / T/Vc _ 3 \ 

16 7 T 2 e V 2 2 N C ) ' 


(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 


An immediate check of our results is the extraction of the correct anomalous dimension, 7 cm, already known in 
the literature for the operator Ocm (Eq. (i37l) ~). with a quark mass and a coupling constant in its definition [2(|. The 
following relation holds between z 1 ’ and 7 cm'- 


~ r, 2 / Dfl.MS . 

7 CM = -2 eg (z 1 ’ + 





(73) 


ryDR, MS -1 , 2 DR,MS 

Z m - 1 +9 z m 


0(g 4 


MQ 1 1 


,DR ,MS _ 


16 7T 2 


(-3 C F ) 


(74) 


3.2. Lattice regularization - MS renormalization 

The computation of the 2 -pt and 3-pt bare Green’s functions of Ocm on the lattice are the most demanding part 
of the present work. This is particularly true for the 3-pt function, since it had to be calculated for arbitrary values 
of the external momenta, qt , of the quarks and gluon. The algebraic expressions involved were split into two parts: 
a) Terms that can be evaluated in the a —> 0 limit: Such terms exhibit a very complicated dependence on qi , even for 
zero quark masses, involving Spence functions. These functions constitute a part of the regularization independent 
renormalized Green’s functions, b) All remaining terms: These are divergent as a —> 0, however their dependence on 
qi,m is necessarily polynomial. Our computations were performed in a covariant gauge, with arbitrary value of the 
gauge parameter a. Given that some of the operators which mix with Ocm contain powers of the quark masses, we 
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have kept these masses different from zero throughout most of the computation; it is only in the final expressions for 
Zi that we set m —> 0. 

For the algebraic operations involved in evaluating Feynman diagrams, we make use of our symbolic package in 
Mathematica. A brief description of the computation of a Feynman diagram can be found, e.g., in Ref. [2l| and 
references therein. The algebraic expressions for each Feynman diagram typically involve ~ 10 5 terms at intermediate 
stages. The requirements in terms of CPU time, both for algebraic manipulation and for numerical integration of 
momentum loop integrals, were rather modest as compared to human effort: a total of ~ 4 months on a single core 
CPU was required. 

The computation on the lattice is performed in the twisted basis (ip', ip'), and thus, before comparing with the 
results in DR, we must rotate to the physical basis (ip, ip). This rotation amounts to the following transformation of 
the fermion field: 


ip' = e -i ip , (75) 

ip' = e-^^ip. (76) 

The rotation of the 2-pt Green’s function is therefore: 

(Ip o iP) amp = e~* 4 (ip' O ip') amp e- 1 4 * , (77) 

and similarly for the 3-pt Green’s function. 

We will make use, once again, of Eqs. ([35]) - (l36l) . with MS being the renormalization scheme; however, the 
regularization will now be the lattice. The above equations now take the form: 

E ((^- MS )- i ) 1 .(v»o i ?)L P (78) 


and 


(iPOuPA^X = ^ S (^ MS ) 1/2 E((^ MS )- 1 ) li <^^^ mp . (79) 

i=l 

The left-hand sides of the above equations are known from the calculations in DR, see Subsection 13.11 The bare 
lattice Green’s functions in these equations contain terms which diverge in the limit a —> 0; these divergent terms 
have a form similar to Eqs. (l38l) and (l39l) . with two differences: 

• 7 - log (a 2 ) 

contributions: 

Ps (r d 75 <h + r s #2 75 ) + Ps i (rd m d + r s m s ) 75 + p w ■ 1 (80) 

Rsg(r d - rsj'ish, • (81) 

These contributions lead to mixing with the lower dimension operators On, 0 12 and O 13 defined in Eqs. HMD - 

m- _ _ 

The renormalization functions Z^’ MS (Z^ MS ) for the quark (gluon) field, as well as Zpj ,MS , Zpp ,MS , were only 
partially available in the literature; we computed them for a general covariant gauge, using the Symanzik improved 
gauge action for different sets of values for the Symanzik coefficients. These results are presented in Appendix [B] in 
the RT renormalization scheme along with conversion factors to the MS scheme. 

Renormalizability of the theory implies that the difference between the one-loop renormalized and bare Green’s 
functions must only consist of expressions which are polynomial in qi, m ; in this way, the right-hand sides of Eqs. CB 
- (El can be rendered equal to the corresponding left-hand sides, by an appropriate definition of the (<&- and m- 

independent) renormalization functions Z^’ MS . These differences can be written as follows: 

(ip O! VS - (ip O ! p = 9 2 (z^ - (ip 0 1 ip )tree - ^ ^^ ^ ^tree (82) 


• There are additional O 


, o 


in (iP O i ip) 


amp 


in {ipOiipA v )^ mp lPI : 


13 


and 




s 2 (+ MS + i+ MS + - + MS )» o. 

13 

- Y,Zi MS ^O^A v ) tTee . 

i =2 


(83) 


Indeed, we have checked explicitly the polynomial character of Eqs. (1521) - (1551) . This check is quite nontrivial, 
especially for Eq. (1831) . since both the bare and renormalized Green’s functions, taken individually, exhibit a very 
complex dependence on the momenta qi. The left-hand sides of Eqs. (1521) - (1551) have the same tensorial form as 
Eqs. (1551) - (1551) , respectively, but with the additional contributions of Eqs. (1501) - (1511) . 

Each tensorial structure (multiplying p\ — pio, R\ — R 5 ) will provide an equation; the set of these equations (a 
total of 15) can be solved for the 13 mixing coefficients Z t . Two of the equations serve as consistency checks and the 
remaining 13 lead to a well determined system. Upon solving all equations we obtain for the Iwasaki gluon action 
(see Appendix [A] for other gluon actions we have considered) 


z i 


ryL ,MS 
Z 2 


ryL, MS 
Z 3 

ryL, MS 
Z 4 


ryL, MS 
Z 5 


ryL, MS 
Z 6 


ryL, MS 

Z 7 


ryL,MS 
Z 8 


ryL, MS 
Z 9 


ryL, MS 
^10 


ryL,MS 
Z 11 


ryL,MS 
^12 


ryL,MS 
^13 


16 7 r 2 


N c ( -7.9438 + - log (a 2 p 2 ) j + — ( 4.4851 - - log (a 2 p 2 ) 


9 2 C f 
16 7 r 2 

0 

0 


[4.5370 + 6 log (a 2 p 2 )] 


9 


16 7 r 2 

0 


N c (4.2758 - | log (a 2 p 2 )) + (-3.7777 + 3 log (a 2 /2 2 )) 


7 l, ms 


9 2 C f 
16 7 r 2 


(—3.7760) 




I.,MS 


g 2 <?F 

16 7T 2 

1 g 2 C, 


a 16 7T 2 


[3.7777-3 log (a 2 /2 2 )] 


(-3.2020) 


-Z- 


L, MS 
11 

1 .9 2 Of 
a 2 16 7r 2 


(36.0613) . 


(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 


In these equations, p is the MS renormalization scale which appears in {ipO\ and (ifiO 1 V’-'^amp by virtue of: 

g = p~* ( Z g DR } g = p( 47 r/f°y/ 2 . 

The above results for Z^' MS - zfy MS are independent of the choices r s = ±1, = ±1. There is also a small system¬ 

atic error originating from the numerical estimation of lattice integrals, however it is much smaller than the displayed 
accuracy of the results. It is important to emphasize that the coefficients Z\\—Z 13 , which control the mixing with lower 
dimension operators, may receive also non-perturbative contributions proportional to [( 1 /a) exp (—1 /( 2 /? 0 < 7 o))]* ~ ^ 
(with k = 1,2 for Z 1112 and Z 13 respectively) [22| . For this reason, a proper subtraction of the mixing with operators 
On — O 13 must be implemented in a non-perturbative way [ 12 L1131. Il 8 l|. 

If one wants to renormalize in an (appropriately defined) RT scheme, the calculation in DR is not necessary: it 
suffices to compute the bare Green’s functions on the lattice. In this case the left-hand sides of Eqs. (1551) - (1551) . for 
particular values of the external momenta, are dictated by the RT renormalization conditions. 
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The conversion factor between the RI' and the MS scheme will actually be a 
Since this matrix is regularization independent, one may compute it through: 


£)MS = qRI',MS qRI 1 (jRl \MS _ 


'R = ^ R ! 

Thus, in RI', the mixing coefficients read (in matrix notation): 

gL,RI' gL,MS qRI',MS 




(13x13) matrix in this case: CJ* 1 ,MS . 
Z D R,R!'' ( 97 ) 

(98) 


4. NON-PERTURBATIVE DETERMINATION OF THE POWER-DIVERGENT MIXING 

COEFFICIENTS 

In this Section we present the non-perturbative determination of the coefficients Z \ 3 and Z 12 describing the power- 
divergent mixings of the chromomagnetic operator with the scalar and pseudoscalar densities, respectively. We use 
lattice QCD simulations with the gauge configurations produced by ETMC with four ffavors of dynamical quarks 
(Nf = 2 + 1 + 1 ), which include in the sea, besides two light mass degenerate quarks, also the strange and charm 
quarks with masses close to their physical values [HQ- 

Due to the equations of motion some of the operators 0± - 0 13 do not appear in the calculation of on-shcll matrix 
elements. The remaining ones, namely 0 1 , 02 , 03 , 0 4 , O 12 and 0 13 , are present and it is therefore crucial to 
have a reliable estimate of the corresponding mixing coefficients. For operators of the same dimensionality as the 
chromomagnetic one, i.e. 0\ - O 4 , our perturbative one-loop results [see e.g. Eqs. (184H87[) ] are expected to provide 
a satisfactory estimate. However, as already discussed in Section [TJ for the mixing coefficients of the operators with 
lower dimensionality 012 and 0 13 , which are power-divergent, perturbation theory is expected to provide only a 
ballpark estimate 

In order to achieve a non-perturbative estimate of the mixing coefficients of the operators 0 12 and O 13 we impose 
the following renormalization conditions El 

lim (n\Of\K) = 
m s , H) 

m 1 K) ms}Tnd = 

where the pion and kaon states are taken to be at rest 6 . Note that in this Section the operators 0; are the bare 
(local) lattice versions of the operators introduced in Section [2] [see Eqs. (11911311) 1], In Eqs. (|99I100|) we have factorized 
out explicitly the power divergence of the mixings with the operators 0 12,13 and we have introduced the renormaliza¬ 
tion scale independent [H] mixing coefficients 012 , 13 , which can be written as appropriate combinations of operator 
renormalization constants (see Eq. (11051) below). The renormalization conditions (l99l) and (11001) are valid up to cutoff 
effects, which are 0(a 2 ) for Eq. (1M1) and 0(a) for Eq. (HOOD . 

The first condition (1991) requires that the renormalized chromomagnetic operator 0^ has vanishing on-shell matrix 
elements in the SU(3) chiral limit. This is also consistent with ChPT, which for the matrix element of 0f = 167 t 2 Q+ 
predicts at LO [see Eqs. (BED ] 

11 M 2 

(ni(p„)\Of\K>(p K )) = —& --*—(*>„■ p K )B C MO , (101) 

2 m s + TOd 

where = 1, C° = —\/\/2 and Bcmo is a SU(3) ChPT low-energy constant. 

The second condition (HOOD imposes that in the continuum limit the parity violating matrix element (0|Of \K) is 
identically vanishing. 

As usual, in the lattice version of the chromomagnetic operator (fl9l) the gluon tensor G Ml/ is replaced by its clover 
discretization P Ml/ , namely [Hj 


^- lim (Ti-IOr - ^0 13 \K) = 0 , (99) 

Zj\ m s , m d —>0 a 

^(O|0i - Gig - — 0 12 \K) m ^ md = 0 , (100) 

Zi a z a 


0 i = go , 


(102) 


Such a choice is motivated by the fact that, as already observed, the matrix elements of the operator O 4 = □ 


( 4’si’d. ) 


are proportional 


to the squared four-momentum transfer between the external hadronic states. This implies that for external 7 r- and A'-mesons at rest 
(or, more precisely, in the case of mesons with equal spatial momentum) the operator O 4 has vanishing matrix elements in the SU(3) 
chiral limit and therefore it does not affect the determination of the mixing coefficient C 13 from Eq. (55} . The same happens for the 
contributions of the operators O 2 and O 3 , since they are directly proportional to the quark masses. 
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where 


Puv(x) = 


4 

\ ' 


4a 2 “ 2*5o - 


p j(x) - PUx) 


(103) 


and the sum is over the four plaquettes Pj{x) in the fi-v plane stemming from x and taken in the counterclockwise 
sense. 

Our lattice setup is the same as the one adopted in Ref. 0] for the determination of the urn down, strange and 
charm quark masses. The fermions were simulated using the Wilson Twisted Mass Action [|, 0] which, at maximal 
twist, allows for automatic (H(a)-improvement 0,0. In order to avoid the unphysical flavor mixing in the strange 
and charm valence sectors we adop ted the non-unitary set up described in Ref. 0, in which the valence quarks 
are regularized as OS fermions |27| . For the gauge links we simulated the Iwasaki action 0], because it proved to 
facilitate simulations with light quark masses allowing to bring the simulated pion mass down to approximately 210 
MeV 00| . 

The details of the ETMC gauge ensembles are collected in Table IIII1 where the number of the gauge configurations 
analized (N c f g ) corresponds to a separation of 20 trajectories. At each lattice spacing, different values of the light 
sea quark mass were considered. The up and down quark masses were always taken to be degenerate and equal 
in the sea and valence sectors ( m u = md = rnt). The masses of both the strange and the charm sea quarks 
were fixed, at each /3, to values close to the physical ones 0- We simulated quark masses for the light sector in 
the range 3 m yhys < mt < 12 while we used three values of the valence strange quark mass in the range 

0.7 m p s hvs < m s < 1.2 m phys in order to interpolate to the physical strange quark mass 0. The simulated pion 
masses cover the range ~ 210 4- 450 MeV. 


ensemble 

P 

V/a 4 

ttf^sea — 

a/J, a 

a/* <5 

Ncfg 

afis 

A30.32 

1.90 

32 3 x 64 

0.0030 

0.15 

0.19 

150 

0.0145,0.0185,0.0225 

A40.32 



0.0040 



100 


A50.32 



0.0050 



150 


A60.24 

1.90 

24 3 x 48 

0.0060 

0.15 

0.19 

150 


A80.24 



0.0080 



150 


A100.24 



0.0100 



150 


1325.32 

1.95 

32 3 x 64 

0.0025 

0.135 

0.170 

150 

0.0141,0.0180,0.0219 

1335.32 



0.0035 



150 


1355.32 



0.0055 



150 


1375.32 



0.0075 



80 


1385.24 

1.95 

24 3 x 48 

0.0085 

0.135 

0.170 

150 


1315.48 

2.10 

48 3 x 96 

0.0015 

0.12 

0.1385 

60 

0.0118,0.0151,0.0184 

1320.48 



0.0020 



100 


1330.48 



0.0030 



100 



TABLE III: Details of the gauge ensembles and values of the simulated sea and valence quark bare masses used in this work 
(after Ref. M)- 

Quark propagators with different valence masses are obtained using the multiple mass solver method j29l . [30| , which 
allows to invert the Dirac operator for several quark masses at a relatively low computational cost. The statistical 
accuracy of the meson correlators is significantly improved by using the “one-end” stochastic method [0, which 
includes spatial stochastic sources at a single time slice chosen randomly. Statistical errors are evaluated using the 
jackknife procedure. 

In our lattice setup both the scalar and the pseudoscalar non-singlet densities renormalize multiplicatively and, 
choosing in particular r s =rd—r (see below), we have 

Ol 3 = ZpOu = Zp ^ s ip d , 0? 2 = 1^012 = + HdW.T&'I’d , (104) 

Zj p Zj p 

where fd = Zp ■ m is the (twisted) bare quark mass, while Zp and Zs are the renormalization constants of the 
(pseudo)scalar densities computed using the RI'-MOM scheme in Ref. 0]. Therefore, the mixing coefficients C 12 and 
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C 13 introduced in Eqs. (I99H100D are related to the operator renormalization constants via 


C13 


= Z 13 Zp 


C12 _ „ Zs 

- = ^ 12 ^— 

CL Zp 


( 105 ) 


In order to extract the coefficient C 13 from the condition (1991) we have computed for each gauge ensemble the 3-point 
meson correlators defined as 


CiM = (0\P 5 n (y)O J (x)P 5 K Hz)\0)S t , {tx ^ ) Sp Aty _ tz 


L e 


(106) 


x,y,z 


where j = 1 or j = 13, P 3 (x) = iip d (x)'Ys’ipui.x) and P^{x) = iil> s (x) r y 3 il> u (x). Our choice r u = — r d = — r s guarantees 
that the two valence quarks in the pion and kaon mesons have always opposite values of the Wilson parameter. In 
this way the pion and kaon states in the matrix elements of the renormalized chromomagnetic operator ( 7 r|(H^|Ar) 
have good scaling and chiral properties. 

At large euclidean time the correlation function is dominated by the contribution of the lightest states and one gets 


C J 3 (t,t') 


Z- n 


Z* K 


t^a, (t'-t)>Q 2 M n 2 M k 


{AOj\K) 




(107) 


with Z n[K) = (0 \P^ k \o)\tt(K)). 

From the large time behavior of the 3-point correlators corresponding to the chromomagnetic and scalar density 
insertions one can construct the following ratio 


Cl (MO (MOi\K) 

C 3 13 (M0 t » n , <7T|M 13 |^> 


mg), 


(108) 


where the first two quark masses are those involved in the transition and the third one is the spectator valence quark 
mass, which is taken to be equal to the light sea quark mass. According to Eq. (l99l) . the ratio R\ 3 {m s ,mi;mi) 
provides in the SU(3) chiral limit an estimate of the mixing coefficient C 13 at each value of the lattice spacing, namely 

ci 3 = lim a 2 Ri 3 (m s ,mf,me) . (109) 

m s , m^—^0 


Note that the dimensionless quantity a 2 Ri 3 (m s , mi; mi) is extracted directly from the ratio (11081) computed in lattice 
units. 

In order to perform the chiral limit we start by computing the ratio C 3 (t,t' = T/2)/C 3 3 (t,t' = T/ 2) in the 
degenerate case m 3 = mi for all the gauge ensembles of Tabic M- which will be referred to as the tt —> n channel. 
In this channel the mixing with the operator O 4 is absent and the mixing with the operator O 12 is linear in the light 
quark mass, while the mixings with the operators O 2 and 0 3 are quadratic in the quark mass. 

The results obtained for the ensembles corresponding to the lightest pion mass at each of the three lattice spacings 
are presented in Fig. [5] It can be seen that the ratio C 3 (t,T /2) /C 3 3 (t,T /2) exhibits nice plateaux. In what follows 
we consider two different choices for the time intervals adopted to extract the values of a 2 R 33 (mi, mi; mi) from the 
ratio C 3 (t,T/2)/C 3 3 (t,T/2). The first choice, which will be referred to as the short plateaux (SP), corresponds to the 
time intervals [t m i n ,T/2 — t m i n ], where t m i n is the time distance at which the pion ground state starts to dominate 
the corresponding 2-point correlator according to the analysis carried out in Ref. |25l | (see the horizontal solid lines 
in Fig. [5]). For the second choice, which will be referred to as the long plateaux (LP), the time intervals are extended 
from 6 a to T /2 — 6 a, as shown in Fig. [5] by the horizontal dashed lines. 

The results obtained for a 2 Ri 3 (mi,mi',mi), adopting the LP choice for the plateaux, are collected in Fig. [ 6 ] and 
show a linear dependence on the light quark mass, as expected at leading order in the quark mass expansion. Therefore 
we fit the data at each lattice spacing using a linear ansatz of the form a 2 R± 3 (mi , mi) = C 13 + A ■ am, where the 
parameters C 13 and A are determined by minimizing the y 2 variable. Similar results (with larger statistical errors) 
hold as well for the values of a 2 Ri 3 (mi, mi; mi) obtained adopting the SP choice for the plateaux. 

As a check of the uncertainty related to the chiral extrapolation, we have also computed the ratio a 2 R\ 3 {m s , m s ; mi) 
using for m s the values corresponding to the valence strange (bare) quark mass reported in Table Hill With respect to 
Eqs. (I108I) - (I109D we replace the d quark in the transition by a strange-like quark s' with mass m s ' = m s and ry = r s . 
We refer to this channel as the K —> K one. The results for the ratio a 2 R\ 3 {m s ,m s ;mi) obtained for the ensembles 


7 


Precisely, we replace in Eqs. 1 X 021 ) -CM) the strange quark s by a light quark d! with mass = vni and r^r = In this way we 
guarantee the absence of disconnected contributions and keep the discretization effects in the squared pion mass of order 0 (a?m#). 
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FIG. 5: The ratio of 3-point correlators C\(t,t') / C\ 3 in lattice units for the gauge ensembles A40.24 (a), B25.32 (b) and 
D15.48 (c) versus the insertion time t for a fixed time separation t! = T/2 between the source and the sink. The horizontal 
solid (dashed) lines correspond to the central values and the errors obtained adopting the SP (LP) choice of the plateaux (see 
text). 


at f3 = 2.10 are reported in Fig. [3(a), where it can be seen that the data at fixed value of the light quark mass 
can be fitted adopting a quadratic ansatz in the strange quark mass. In Fig. (3b) the SU(3) chiral point is finally 
reached by performing a linear fit in the light quark mass and the result is compared with the one corresponding 
to the 7 r — > n channel [see Fig. [ 6 ] at /3 = 2.10]. A good agreement between the two channels is obtained within one 
standard deviation. Similar results are obtained at (3 = 1.90 and 1.95. 

Our determinations of the mixing coefficient C 13 obtained from the n —> 7 r and K —> K matrix elements, adopting the 
two choices SP and LP for the plateaux, are presented in Table ITVl A very good agreement is found between the results 
obtained in the two channels, while small differences (within at most ft! 1.5 standard deviations) are visible between 
the values corresponding to the SP and LP choices for the plateaux. This comparison will be useful to quantify 
the uncertainties due to the subtraction of the mixing with the scalar density in the study of the (renormalized) 
chromomagnetic operator matrix elements [l7j, I3- Note that the uncertainties on C 13 are of the order of 0.01% in 
the case of the LP choice, while they do not exceed the level of 0.1% in the case of the SP choice. 

We now want to compare our non-perturbative results of Table hy i with the predictions of perturbation theory at 
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FIG. 6: The ratio a 2 Ri 3 {mt,me\mt), extracted using the LP choice for the plateaux, versus the (twisted) quark bare mass 
afu = Zpame for each value of the lattice spacing. The solid lines are the results of linear fits in age applied to all data. The 
values of the \ 2 variable (divided by the number of degrees of freedom) are reported in each inset. The diamonds represent 
the values of the mixing coefficient C13 obtained in the chiral limit [see Eq. (11091) 1. 


p 

7r —► 7r channel 

K —> K channel 


SP 

LP 

SP 

LP 

1.90 

0.89769 (17) 

0.89710 (11) 

0.89752 (24) 

0.89716 (10) 

1.95 

0.87687 (36) 

0.87627 (12) 

0.87687 (38) 

0.87632 (13) 

2.10 

0.81646 (78) 

0.81675 (08) 

0.81635 (61) 

0.81677 (08) 


TABLE IV: Values of the mixing coefficient C13 obtained as the chiral limit of the data on the ratios a 2 Riz{rnt, me', me) for the 
7r —> 7r channel and a 2 R\ 3 (m s , m s ; me) for the K —► K channel, extracted at three values of the lattice spacing using the SP 
and LP choices for the plateaux. 


one loop obtained in the previous Section. Since the mixing coefficient Z 13 starts already at order 0(g 2 ), while both 
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FIG. 7: (a) The ratio a 2 Ri 3 (m a , m s \me) in case of the K —► K channel, adopting the LP choice for the plateaux, versus 
the strange quark bare mass ag a = Zpam a for fixed values of the light quark mass age = Zpami corresponding to the three 
ensembles B15.48, B20.48 and B30.48 at = 2.10. The dashed lines are quadratic fits in ag s applied to the data, (b) The 
ratio a 2 f?i3(0, 0; me) for the K —>• K channel versus the light-quark bare mass age. The solid line is a linear fit in age and the 
empty square is the corresponding value of C13. The empty circle is the result for C13 obtained from the 7r —» 7 r channel [see 
Fig. [f3 at /3 = 2.10]. 


Z\ and Zp start at order 0(1), the one-loop perturbative term for C 13 , defined in Eq. (I105D . is the same as the one 
of a 2 Z 13 , namely in the case of the ETMC action (see Eq. (i96l) f 


c i3 looP PT = 36.0613 . (110) 

In Fig. [H using g 2 = g% = 6//3, the non-perturbative results for ci 3 ( 167 t 2 )/(< 7 2 Cf), obtained at three values of 
the lattice spacing using the LP choice for the plateaux, are compared with the corresponding perturbative result 
from Eq. (I110D . It can be seen that the non-perturbative determinations of C 13 differ less than ~ 10% from the 
perturbative predictions at one loop at all values of the lattice spacing. We expect however that in order to get a 
reliable determination of the renormalized CMO matrix elements a high-precision determination of C13, at the level 
of 0 . 1 % or better, will be required ehi. 

Let us now discuss the non-perturbative determination of the mixing coefficient C12 of the chromomagnetic operator 
with the pseudoscalar density. To this end we make use of Eq. (11001) and of our accurate non-perturbative results for 
the mixing coefficient C 13 . We anticipate that the O(a) terms affecting the r.h.s. of Eq. (11001) . which are unavoidably 
relevant in the parity violating matrix elements {0\O\\K) and (O|0i3|.A}, turn out to be numerically competitive for 
our values of the lattice spacing with the contribution of the power-divergent mixing with the pseudoscalar density 
012 • This finding can be expected also from the smallness of the one-loop perturbative estimate of C12 (see Eq. (ED) 
with respect to the corresponding result (l96l) for the mixing coefficient C13. 

We start by computing for each gauge ensemble the 2-point meson correlators defined as 

= ^3 E (°l P 5 K (x)P 5 K \z) |0> 5 tXtx - tz) , 

x,z 

clit) = ^ (0| Oj{x)P^\z) |0) j 

x,z 

where j = 1 or j = 13 and P^{x) = i'il> s (x)'y 5 il’d(x ) with r s = rj (see Eq. (11041) 1. At large time distances one has 


( 111 ) 

( 112 ) 


Cf(i) 


t»a 


cm 


t»a 


\Z K \ 2 M K t 

2 M k 

Z K (Q\Oj\K) Mr<t 
2 M k 


(113) 
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FIG. 8: The quantity Ci3(167r 2 )/(<? 2 Cf) versus the coupling g 2 /(An) = <7o/(4tt) = 3/(2n/3) calculated non-perturbatively at 
three values of the lattice spacings using the LP choice for the plateaux in the n —>• n channel (see Table nvj The horizontal 
dotted line is the prediction of lattice perturbation theory at one loop corresponding to Eq. cm The solid line corresponds 
to a linear fit in g 2 /(An) and the full circle is the corresponding extrapolated value at g 2 = 0. 


with Zk = (0|Pg f |ii'}. We stress that, since r s = r f j, the two valence quarks in the kaon have the same value of 
the Wilson parameter and therefore the squared meson mass M^ differs from its continuum counterpart by terms of 
order 0 (a 2 ), which do not vanish in the chiral limit. 

Then, using our non-perturbative results for C 13 , the one-loop perturbative estimates for Z\ (see Eq. l|84Jl ) and the 
non-perturbative renormalization constant Zg from Ref. [25| , we have computed the following ratio 


Pi 2 (m^,m s ;a) 


1 

ZMQIP^K) 


(0|Oi| K) 


-^(0|Oi 3 | K) 


(114) 


both in the degenerate case m s = me and for the three simulated values of m s reported in Table IIIII For the 
dimensionless quantity aroRi 2 (m s ,me', a), where is the Sommer parameter, the continuum limit expectation is a 
simple linear dependence on the sum of the light and strange (renormalized) quark masses, viz. 


ar 0 R 12 (me,m s ;a) = c 12 


Z P 

ZsZ\ 


ro(mg + m s ) + 0(a 2 ) , 


(115) 


where the factor Zp / (ZgZ\) is almost constant for our three /3-values, namely Zp / (ZsZ\) = {1.285,1.275,1.260} for 
/3 = {1.90,1.95,2.10}. Note that the ratio Zp/Zs is both scheme and renormalization scale independent, while the 
renormalization constant Z\ carries the scheme and renormalization scale dependence of the l.li.s. of Eq. (11151) . Thus 
the coefficient C 12 is both scheme and renormalization scale independent (see Ref. 0). 

However, as anticipated, the contribution of the lattice artefacts in Eq. (11151) is not negligible. Beyond trivial terms 
proportional to a 2 and a 2 (me + m s ), coming from the mixing of Oi with dimension-6 parity-odd operators, one should 
also consider that in the twisted-mass approach the Symanzik expansion of correlators like C 2 (t) (see Eq. (11121) 1 may 
contain chirally enhanced cutoff effects that are described by the (space-time integrated) insertion of the parity-odd 
local operator £ 0 dd = £5 + 0(a 3 ), connecting the vacuum with the one-pion state [32;, 33}. In particular the matrix 
element of the chromomagnetic operator (0|Of \K) may receive a cutoff effect proportional to (0|£5|7r)(7r|C>f \K)/M 2 . 
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Therefore, taking into account that at maximal twist we have {0 |jCs |7t) oc ame [321 ] • the following fit has been performed 


ar 0 R 12 (m s ,me;a) 


Zp \ 

Cl2 inr + ^ 12 ° 2 ) r °( m ^ + TOfi ) +ho a 2 + h' 0 a 4 

2 Mk u 2 (Mk ~ 
h c a r 0 + h s a r 0 m e -- , 


(116) 


where C 12 , c?i 2 , ho, h r 0 , he and h$ are fitting parameters and the last two terms come from the expected mass 
dependence of the matrix elements of the chromomagnetic operator, (Tt\Oi\K) oc M^Mk, and of the operator 
Oa = □ (ips'i/Jd), (7t|D ('tp s ipd) | K) oc (. M k - Mp 1 . 

The kaon mass Mk appearing in Eq. (11161) is directly extracted from the lattice correlator (t) (see Eq. (11111) ) and 
contains 0(a 2 ) terms, which make it heavier than its continuum counterpart. Instead the pion mass appearing in 
Eq. (11161) may receive a tower of cutoff effects by multiple insertions of the parity-odd local operator £ 0 dd> so that it 
is not known a priori whether it may be lighter or heavier than its continuum counterpart. Therefore we have adopted 
for the pion mass two choices differing by 0(a 2 ) terms, namely the OS pion M° s [27]], which coincides with Mk in 
the mass-degenerate case m s = mi, and the twisted-mass neutral pion M^ 0 M , which for our lattice setup turns out 
be lighter than the charged one. For the values of Mj 0 M we have used the lattice results reported in Ref. j34|. 

We have applied the fitting function (11161) to the lattice data for ar 0 Ri 2 (m s ,mf, a) and found that the use of M t 0 m 
produces smaller values of x 2 with respect to the OS choice M° s . Moreover we have tried also to insert an extra 
cutoff effect in the pion mass of the form M 2 = [Mj M ] 2 + ca 2 . However, the resulting value of the fitting parameter 
c turned out to be largely compatible with c = 0 and no improvement in the \ 2 value was found. 

The good quality of our best fit (corresponding to % 2 /d.o.f. ~ 0.6) is shown in Fig. [9] In particular it can be seen 
that Eq. (11161) is able to take properly into account the dependence on both the lattice spacing (see panel (a)) and 
the light-quark mass (see panel (b)). 



FIG. 9: The dimensionless quantity aroRi 2 (me, m s ; a) versus the (renormalized) quark mass (mi + m s ), calculated using the 
values of ro/a from Ref. [25[ and the non-perturbative result for C13 corresponding to the n —> n channel and to the LP 
choice for the plateaux (see Table HVll . In (a) the dots, squares and diamonds correspond to the results obtained for the gauge 
ensembles A30.32, B25.32 and D20.48, respectively, which share a (renormalized) light-quark mass me ~ 12 MeV. In (b) the 
results corresponding to the ETMC ensembles with /3 = 1.95 (fixed value of the lattice spacing) and various values of the 
light-quark mass me (see the inset) are presented. The dashed lines represent the best fit curves according to Eq. (11161) . 


The final non-perturbative result for the mixing coefficient C 12 is 

Ci 2 = 0.035 (20) , 


(117) 


which is approximately 40% of the value expected from one-loop perturbation theory for the ETMC action, namely 
(see Eq. (I951D 

c i-ioo P pt _ 2__£ 3.202Q = { 0 . 08 54, 0. 0832 , 0.0772} for /? = {1.90,1.95,2.10} . (118) 
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The uncertainty in the non-perturbative determination (11171) of c 12 is « 60%, and this justifies a posteriori that 
the g 2 -dependence of the mixing coefficient C 12 has been neglected in the present analysis. 

We note, in conclusion, that the smallness of both the perturbative (I118D and non-perturbative (11171) determinations 
for C 12 implies that in our lattice formulation the subtraction of the power-divergent mixing with the pseudoscalar 
density is not going to play a crucial role for the numerical determination of the CMO matrix elements. 


5. SUMMARY 

The study of the chromomagnetic operator on the lattice has been hampered up to now by the complicated pattern 
of operator mixing. We identified these operators based on the symmetries of the regularized theory, which has been 
conveniently chosen to be twisted-mass LQCD (at maximal twist) with the Iwasaki gluon action. 

There are mixings with lower dimensional operators (which are power divergent), as well as with gauge non-invariant 
operators. We have computed all relevant mixing coefficients to one loop in lattice perturbation theory; this has 
required the calculation of both 2-point (quark-quark) and 3-point (quark-quark-gluon) Green’s functions at nonzero 
quark masses. We have calculated all the elements of the mixing matrix that is relevant for the renormalization of 
the chromomagnetic operator at one loop in lattice perturbation theory. 

For the first time the l/a 2 -divergent mixing of the chromomagnetic operator with the scalar density has been 
determined non-perturbatively with high precision (see Table El). The 1/a-divergent mixing with the pseudoscalar 
density, which is peculiar of the twisted-mass formulation, has been also calculated non-perturbatively (see Eq. (11171) ) 
and found to be smaller than its one-loop perturbative estimate (I118D . We have carried out the QCD simulations on 
the lattice using the gauge configurations produced by ETMC with Nf = 2 + 1 + 1 dynamical quarks, which include 
in the sea, besides two light mass degenerate quarks, also the strange and charm quarks with masses close to their 
physical values. Three values of the lattice spacing between ~ 0.6 and ~ 0.9 fm and pion masses in the range 210-1-450 
MeV have been considered. 

The results presented in this paper, which determine the mixing pattern of the chromomagnetic operator, are an 
essential ingredient for the determination of the (renormalized) CMO matrix element between pion and kaon states, 
whose calculation is in progress. Preliminary results have been presented in Ref. na and the final ones will be the 
subject of a forthcoming publication [l8[ . 
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Appendix A: Mixing coefficients Z, 

In this Appendix we present our results for the mixing coefficients, Zi (i = 1,...,13) in the MS scheme, for 
the following gluon actions: Wilson, tree-level Symanzik, Tadpole Improved Liischer-Weisz (TILW, at /3 Co = 8.30; 
/3 = 2 N c /g 2 ), Iwasaki and Doubly Blocked Wilson (DBW2). The values of the Symanzik coefficients corresponding 
to these actions are collected in Tabic IVl 


Coefficient 

Wilson 

Tree-level Symanzik 

TILW 03 co = 8.30) 

Iwasaki 

DBW2 

co 

1 

5/3 

2.386978 

3.648 

12.2688 

Cl 

0 

-1/12 

-0.159128 

-0.331 

-1.4086 

C 2 

0 

0 

0 

0 

0 

C 3 

0 

0 

-0.014244 

0 

0 


TABLE V: Symanzik coefficients for various choices of gluon actions. 
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Our calculation has been performed in an arbitrary covariant gauge. All the mixing coefficients Z* (* = 1,..., 13) 
in the MS scheme are gauge independent. The generic forms of the mixing coefficients are 


7 l ,ms _ 


Z. 


z, 


z 


= 1 + 


16 7r 2 


No ( ei,i + - log (a 2 ^ 2 ) ) + — ( ei, 2 - - log(a 2 ^ 2 ) 


,ms g 2 Cf 

16 7T 2 

0, 

0, 


1., MS 
2 

■L,MS 

3 

•L,MS 

4 

1., MS 

5 

L,MS 
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16 7T 2 

0, 


[e 2 + 6 log (a 2 At 2 )] , 


N c ( e 5 ,i - - log(a 2 /i 2 ) ) + — (e 5>2 + 3 log(a 2 /i 2 )) 


_ yL, MS 


yL, MS 

z 7 


L ,MS _ g 2 Cf 
16 7T 2 
V L,MS 

L, MS _ ^5 

9 — o 

L.MS 

10 


( e 8) , 


L,MS 

11 

L, MS 
12 

I., MS 


g 2 C F 
16 7T 2 

i g 2 c F 

a 16 7T 2 


[—e 5 ,2 - 3 log (a 2 /r 2 )] , 
(en) , 


V L, MS 

^11 7 

1 9 2 C f 


(ei3) 


13 a 2 16 7r 2 

where the values of ei, are shown explicitly in Table ED 


(Al) 

(A2) 

(A3) 

(A4) 

(AS) 

(A6) 

(A7) 

(AS) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 


Coefficient 

Wilson 

Tree-level Symanzik 

TILW (Pc o = 8.30) 

Iwasaki 

DBW2 

ei,i 

-16.8770 

-12.8455 

-10.4920 

-7.9438 

-3.2465 

ei ,2 

13.4540 

9.3779 

7.0022 

4.4851 

-0.5102 

e2 

1.9290 

2.7677 

3.4589 

4.5370 

8.5250 

es.i 

5.9806 

5.3894 

4.9311 

4.2758 

2.2834 

e5,2 

-6.4047 

-5.5061 

-4.8014 

-3.7777 

-0.5292 

68 

-4.0626 

-3.9654 

-3.8894 

-3.7760 

-3.4713 

en 

-4.4977 

-4.0309 

-3.6792 

-3.2020 

-1.9216 

ei3 

54.9325 

47.7929 

42.6253 

36.0613 

19.9812 


TABLE VI: Results for the mixing coefficients at one-loop using the MS scheme on the lattice. The finite parts d and e % ,j are 
given for five actions: Wilson, tree-level Symanzik, TILW (/3co = 8.30), Iwasaki and DBW2. 


Appendix B: Perturbative one-loop Renormalization of Z c , , Zm , Za , Z g on the Lattice 

In this Appendix we provide the results of our one-loop calculation for the renormalization functions of the ghost 
field (Z c ), quark field (Z^), gluon field (Za), coupling constant (Z g ), quark mass (Z m ). These functions enter 
the renormalization of the chromomagnetic operator through Eqs. (IB16I) . (1551) . (13S1) . (1551) . (1371) . The computation 
was performed using twisted mass fermions, Symanzik improved gluon action and a general covariant gauge. Here 
we present the results for Wilson, tree-level Symanzik, TILW ((3 cq = 8.30), Iwasaki and DBW2 actions. For the 
extraction of the renormalization functions, we applied the RT scheme at a scale /2. Once we have computed the 
renormalization functions in the RT scheme we can construct their MS counterparts using conversion factors which 
are known (see, e.g., Hi), up to the required perturbative order. 
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The aforementioned renormalization functions are defined as follows 


do 

— Zg <7, 

(Bl) 

c 

= y/Z~ c c R , 

(B2) 

$ 

= 

(B3) 


k, a. 

II 

(B4) 

a 

= a*, 

(B5) 

m 

= Z m m R . 

(B6) 


In the above, Z g actually stands for Z g ’ RI ; similarly for all other Z’s. The renormalization function Z a for the gauge 
parameter receives no one-loop contribution. 


1. Ghost Field Renormalization Z c 


The ghost field renormalization enters the evaluation of Z g (see subsection IBIfl) and it can be extracted from the 


RI' condition 


lim 

a —>-0 




= 1 , 


(B7) 




where Z,^{q, a) is the ghost self energy up to one-loop computed from the diagrams in Fig. [Tol namely 

E R (q,a) = q 2 +0{g 2 ). (B8) 

The generic form of Z R,K1 is 


Z L,m' = j + 9 2 N c r - l 


16vr 2 [ 6c “ L2029a - -j; (3 - a) log (a 2 g 2 ) 

The numerical values of the coefficient e c are listed in Table IVlD for all gluon actions we have considered. 


(B9) 



-►- 



FIG. 10: One-loop Feynman diagrams contributing to the renormalization of the ghost field. A wavy (dotted) line represents 
gluons (ghosts). 


2. Renormalization of Fermion Field (Z^,) and Mass (Z m ) 


In order to obtain the renormalization functions of fermionic operators we also compute the quark field renormal¬ 
ization, Z^, as a prerequisite. 

Zif, is extracted from the RI' renormalization condition on the fermion self energy E ^{q, a) = itf + m + 0{g 2 ), namely 


lim 

a— 


ryL, RI' 


{an) tr (E^(g, a) <j) / {Ai q 2 ) 


= 1. 


(BIO) 


The trace here is over Dirac indices; a Kronecker delta in color and in flavor indices has been factored out of the 
definition of E,('. The Feynman diagrams contributing to E^ are shown in Fig. fi ll Our result for Z^ is 


/7-L.RI' 


1 + 


9 2 C f 

167t 2 


e ^ — 4.7920a + a log (a 2 g 2 ) 


(BH) 


The part of E^ proportional to the unit matrix in Dirac space leads to the value of Z m . Our result for Z m is 


Zt m ' = 1 


9 2 C f 


e m + a - 3 log ( a 2 g 2 ) 


167r 2 L 

The numerical values of the coefficients and e m are listed in Table IVIII 


(B12) 
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FIG. 11: One-loop Feynman diagrams contributing to the renormalization of the fermion field. A wavy (solid) line represents 
gluons (fermions). 


3. Gluon Field Renormalization Za 


The renormalization for the gluon field, Za, can be evaluated from the gluon propagator G^ u (q,a) with radiative 
corrections, namely 


G%(q, a) 


[ <W ~ qyWg 2 q^/q 2 

f n r(aq) a n L {aq) 


(B13) 


where the one-loop contributions to the transverse (IIt) and longitudinal (IIl) parts of the gluon self-energy, 
II= 1 + 0(g 2 ) are obtained from the diagrams of Fig. [T5] The normalization condition is 


lim 

CL —^0 . 


Z A Rl (a^)U T {aq) 


= 1. 


q z —fi z 


(B14) 


Our result up to one-loop is 

„2 


7 l,ri' 


= 1 + 


9 


167T 2 - 


N c ( e A ,i - 0.8863a + -a 2 ) + —e A , 2 ~ 2.1685 N f + [ -N f - —N c + -aN c ) log (a 2 fl 2 ) , (B15) 


13 


(Nf stands for the number of flavors). The numerical values of the coefficients and e A ,2 are listed in Tabic IVTTl 



WWWV^VWWW' 

FIG. 12: One-loop Feynman diagrams contributing to the renormalization of the gluon field. A wavy (solid, dotted) line 
represents gluons (fermions, ghosts). A solid box denotes a vertex from the measure part of the lattice action. 

One may similarly deduce the value of Z a [see Eq. GSD] from the longitudinal part II L(aq) of the gluon self-energy; 
as mentioned before, Z a receives no one-loop contribution. 


4. Coupling constant renormalization Z g 


Z g can be extracted either from the gluon-quark-quark Green’s function, or equivalently from the gluon-ghost-ghost 
Green’s function G^ gc ; we have chosen to compute the latter 8 . It is customary to renormalize the strong coupling 
constant in the MS scheme also when RI' schemes are adopted for the operators. The same choice is made here. 

The corresponding renormalization condition 9 


lim 

ci —^0 


Z L,RI' (2 £RI')1/2 Z L,RI ' G L Sc{q , a) 


q z =H z 


_ /^finite 

^ Acc 5 


(B16) 


8 We have checked, via explicit computation in the Feynman gauge, that the two determinations lead to identical results. 

9 Eq. JB161) is evaluated at vanishing ghost momentum; q stands for the ghost/gluon momentum. 
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where the expression G^" 1 * 6 is required to be the same as the one stemming from the continuum 


lim 
e->0 L 




G finite 
Acc • 


(B17) 


[In the above equation Z RR ’ Rl is required to eliminate only the pole parts of the left-hand side, without additional 
finite terms; hence, it is trivially equal to Z RR ’ MS .] Thus, G®" 1 * 6 is found to be 


/^finite _ i 

u Acc ~ 1 


r/169 3 

1 2 

1 , 

fp , 2 \ 

\ „ r 5 „ T ■ 

[vT2 + T + 

8“ + 

-a log 

V9 2 ) 



The Feynman diagrams contributing to G ASc are shown in Fig. [13] Our result for Z R ’ m is 


z l,rV = 1 + 


16 7T 2 


e g ,i N c + j^-e gt 2 + 0.5287 N f + ( H Nc - ^N f ) log(a" ) 


(B18) 


(B19) 


The numerical values of the coefficients e 9j i and e g ^ are listed in Table lYEl 



FIG. 13: One-loop Feynman diagrams contributing to G A ^ C . A wavy (dotted) line represents gluons (ghosts). 


Coefficient 

Wilson 

Tree-level Symanzik 

TILW (/3 co = 8.30) 

Iwasaki 

DBW2 

e c 

4.6086 

3.7759 

3.2208 

2.5469 

0.9433 

G-ijj 

16.6444 

13.0233 

10.7153 

8.1166 

2.9154 

G-m 

16.9524 

13.6067 

11.4247 

8.8575 

2.9060 

e A ,i 

22.3157 

10.3088 

2.4199 

-7.2464 

-28.5805 

e A ,2 

-19.7392 

-6.6595 

2.0039 

11.8888 

32.2815 

e 9 ,1 

-13.4192 

-6.5831 

-2.0835 

3.4235 

15.6942 

e S,2 

9.8696 

3.3297 

-1.0019 

-5.9444 

-16.1407 


TABLE VII: The coefficients e c , , e m , eA,i, eA ,2 , e 9 ,i and e 9 ,2 for five actions: Wilson, tree-level Symanzik, TILW (/3 Co = 

8.30), Iwasaki and DBW2. 


5. Conversion to the MS scheme 


Each renormalization function on the lattice, Z LRl , may be expressed as a power series in the coupling constant 
g which is already renormalized in the MS scheme (see Section [1}. 

As already mentioned, our one-loop calculations for Z c , ZZ m , Z A and Z g are performed in a generic gauge with 
parameter a RI . The conversion of a RI to the MS scheme is given by 


„Ri 


gL, MS 
7 L, RI' 


y L, MS 
“ A MS 


7 L. RI' 


J A 


(B20) 


Since (Z R,MS /Z R,RI ') = (z£ R,MS /Z® R ’ RI ') = 1 at three loops HH, it follows 


a RI ' = Za 


L,MS 


MS 


L,RI 


7 a = 


„MS 


C A (g^,a^) 


(B21) 
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Since the ratio of Z’s appearing in Eq. (IB21I) must be regularization independent, it may be calculated more easily 
in DR [35|; to one loop, the conversion factor Ca equals 


Ca(9, a) = ——= = 1 

’ ryDRMS 


9 


36(16tt 2 ) 


[(9a 2 + 18a + 97) N c - 407V/] , 


(B22) 


where hereafter both g and a are expressed in the MS scheme. 

Thus, once we have computed the renormalization functions in the RI' scheme we can construct their MS counter¬ 
parts using conversion factors which, up to the required perturbative order, are given by 


C c (g, a) 
C^{g,a) 
^Vn(.9i £*0 


zb Ri' 


7 l, ms 

Ac 

ryL, RI' 

_ 


ryL, MS 



ryL.MS 

Am 


gnu, ri' 


ryDR, MS 
Ac 

ryDR.Rl' 
ryDR, MS 
yDR, RI' 

Am _ . 


ryDR.MS 

Am 


167T 2 


167T 2 


167T 2 


N c 


Cf a, 


Cp (4 + a). 


(B23) 
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